Single-cell reconstruction and mutation enrichment analysis identifies dysregulated cardiomyocyte and endothelial cells in congenital heart disease.

Congenital heart disease (CHD) is one of the most prevalent neonatal congenital anomalies. To catalog the putative candidate CHD risk genes, we collected 16,349 variants [single-nucleotide variants (SNVs) and Indels] impacting 8,308 genes in 3,166 CHD cases for a comprehensive meta-analysis. Using American College of Medical Genetics (ACMG) guidelines, we excluded the 0.1% of benign/likely benign variants and the resulting dataset consisted of 83% predicted loss of function variants and 17% missense variants. Seventeen percent were de novo variants. A stepwise analysis identified 90 variant-enriched CHD genes, of which six (GPATCH1, NYNRIN, TCLD2, CEP95, MAP3K19, and TTC36) were novel candidate CHD genes. Single-cell transcriptome cluster reconstruction analysis on six CHD tissues and four controls revealed upregulation of the top 10 frequently mutated genes primarily in cardiomyocytes. NOTCH1 (highest number of variants) and MYH6 (highest number of recurrent variants) expression was elevated in endocardial cells and cardiomyocytes, respectively, and 60% of these gene variants were associated with tetralogy of Fallot and coarctation of the aorta, respectively. Pseudobulk analysis using the single-cell transcriptome revealed significant (P < 0.05) upregulation of both NOTCH1 (endocardial cells) and MYH6 (cardiomyocytes) in the control heart data. We observed nine different subpopulations of CHD heart cardiomyocytes of which only four were observed in the control heart. This is the first comprehensive meta-analysis combining genomics and CHD single-cell transcriptomics, identifying the most frequently mutated CHD genes, and demonstrating CHD gene heterogeneity, suggesting that multiple genes contribute to the phenotypic heterogeneity of CHD. Cardiomyocytes and endocardial cells are identified as major CHD-related cell types.NEW & NOTEWORTHY Congential heart disease (CHD) is one of the most prevalent neonatal congenital anomalies. We present a comprehensive analysis combining genomics and CHD single-cell transcriptome. Our study identifies 90 potential candidate CHD risk genes of which 6 are novel. The risk genes have heterogenous expression suggestive of multiple genes contributing to the phenotypic heterogeneity of CHD. Cardiomyocytes and endocardial cells are identified as major CHD-related cell types.

Lack of ethnic diversity in single-cell transcriptomics hinders cell type detection and precision medicine inclusivity.

Perhaps one of the most revolutionary next generation sequencing technologies is single-cell (SC) transcriptomics, which was recognized by Nature in 2013 as the method of the year. SC-technologies delve deep into genomics at the single-cell level, revealing previously restricted, valuable information on the identity of single cells, particularly highlighting their heterogeneity. Understanding the cellular heterogeneity of complex tissue provides insight about the gene expression and regulation across different biological and environmental conditions. This vast heterogeneity of cells and their markers makes identifying populations and sub-clusters especially difficult, even more so in rare cell types limited by the absence of rare sub-population markers. One particularly overlooked challenge is the lack of adequate ethnic representation in single-cell data. As the availability of cell types and their markers grow exponentially through new discoveries, the need to study ethnically driven heterogeneity becomes more feasible, while offering the opportunity to further elaborate ethnicity-related heterogeneity. In this commentary, we will discuss this major single-cell limitation particularly focusing on the repercussions it has on disease research, therapeutic design, and precision medicine.

Overlapping pathogenic de novo CNVs in neurodevelopmental disorders and congenital anomalies impacting constraint genes regulating early development.

Neurodevelopmental disorders (NDDs) and congenital anomalies (CAs) are rare disorders with complex etiology. In this study, we investigated the less understood genomic overlap of copy number variants (CNVs) in two large cohorts of NDD and CA patients to identify de novo CNVs and candidate genes associated with both phenotypes. We analyzed clinical microarray CNV data from 10,620 NDD and 3176 CA cases annotated using Horizon platform of GenomeArc Analytics and applied rigorous downstream analysis to evaluate overlapping genes from NDD and CA CNVs. Out of 13,796 patients, only 195 cases contained 218 validated de novo CNVs. Eighteen percent (31/170) de novo CNVs in NDD cases and 40% (19/48) de novo CNVs in CA cases contained genomic overlaps impacting developmentally constraint genes. Seventy-nine constraint genes (10.1% non-OMIM entries) were found to have significantly enriched genomic overlap within rare de novo pathogenic deletions (P value = 0.01, OR = 1.58) and 45 constraint genes (13.3% non-OMIM entries) within rare de novo pathogenic duplications (P value = 0.01, OR = 1.97). Analysis of spatiotemporal transcriptome demonstrated both pathogenic deletion and duplication genes to be highly expressed during the prenatal stage in human developmental brain (P value = 4.95 X 10-6). From the list of overlapping genes, EHMT1, an interesting known NDD gene encompassed pathogenic deletion CNVs from both NDD and CA patients, whereas FAM189A1, and FSTL5 are new candidate genes from non-OMIM entries. In summary, we have identified constraint overlapping genes from CNVs (including de novo) in NDD and CA patients that have the potential to play a vital role in common disease etiology.

Analyzing single cell transcriptome data from severe COVID-19 patients.

We describe the protocol for identifying COVID-19 severity specific cell types and their regulatory marker genes using single-cell transcriptomics data. We construct COVID-19 comorbid disease-associated gene list using multiple databases and literature resources. Next, we identify specific cell type where comorbid genes are upregulated. We further characterize the identified cell type using gene enrichment analysis. We detect upregulation of marker gene restricted to severe COVID-19 cell type and validate our findings using in silico, in vivo, and in vitro cellular models. For complete details on the use and execution of this protocol, please refer to Nassir et al. (2021b).

Single-cell transcriptome identifies molecular subtype of autism spectrum disorder impacted by de novo loss-of-function variants regulating glial cells.

BACKGROUND: In recent years, several hundred autism spectrum disorder (ASD) implicated genes have been discovered impacting a wide range of molecular pathways. However, the molecular underpinning of ASD, particularly from the point of view of 'brain to behaviour' pathogenic mechanisms, remains largely unknown. METHODS: We undertook a study to investigate patterns of spatiotemporal and cell type expression of ASD-implicated genes by integrating large-scale brain single-cell transcriptomes (> million cells) and de novo loss-of-function (LOF) ASD variants (impacting 852 genes from 40,122 cases). RESULTS: We identified multiple single-cell clusters from three distinct developmental human brain regions (anterior cingulate cortex, middle temporal gyrus and primary visual cortex) that evidenced high evolutionary constraint through enrichment for brain critical exons and high pLI genes. These clusters also showed significant enrichment with ASD loss-of-function variant genes (p < 5.23 × 10-11) that are transcriptionally highly active in prenatal brain regions (visual cortex and dorsolateral prefrontal cortex). Mapping ASD de novo LOF variant genes into large-scale human and mouse brain single-cell transcriptome analysis demonstrate enrichment of such genes into neuronal subtypes and are also enriched for subtype of non-neuronal glial cell types (astrocyte, p < 6.40 × 10-11, oligodendrocyte, p < 1.31 × 10-09). CONCLUSION: Among the ASD genes enriched with pathogenic de novo LOF variants (i.e. KANK1, PLXNB1), a subgroup has restricted transcriptional regulation in non-neuronal cell types that are evolutionarily conserved. This association strongly suggests the involvement of subtype of non-neuronal glial cells in the pathogenesis of ASD and the need to explore other biological pathways for this disorder.

Single-cell transcriptome identifies FCGR3B upregulated subtype of alveolar macrophages in patients with critical COVID-19.

Understanding host cell heterogeneity is critical for unraveling disease mechanism. Utilizing large-scale single-cell transcriptomics, we analyzed multiple tissue specimens from patients with life-threatening COVID-19 pneumonia, compared with healthy controls. We identified a subtype of monocyte-derived alveolar macrophages (MoAMs) where genes associated with severe COVID-19 comorbidities are significantly upregulated in bronchoalveolar lavage fluid of critical cases. FCGR3B consistently demarcated MoAM subset in different samples from severe COVID-19 cohorts and in CCL3L1-upregulated cells from nasopharyngeal swabs. In silico findings were validated by upregulation of FCGR3B in nasopharyngeal swabs of severe ICU COVID-19 cases, particularly in older patients and those with comorbidities. Additional lines of evidence from transcriptomic data and in vivo of severe COVID-19 cases suggest that FCGR3B may identify a specific subtype of MoAM in patients with severe COVID-19 that may present a novel biomarker for screening and prognosis, as well as a potential therapeutic target.

Association of cystic fibrosis transmembrane conductance regulator with epithelial sodium channel subunits carrying Liddle's syndrome mutations.

The association of the cystic fibrosis transmembrane conductance regulator (CFTR) and epithelial sodium channel (ENaC) in the pathophysiology of cystic fibrosis (CF) is controversial. Previously, we demonstrated a close physical association between wild-type (WT) CFTR and WT ENaC. We have also shown that the F508del CFTR fails to associate with ENaC unless the mutant protein is rescued pharmacologically or by low temperature. In this study, we present the evidence for a direct physical association between WT CFTR and ENaC subunits carrying Liddle's syndrome mutations. We show that all three ENaC subunits bearing Liddle's syndrome mutations (both point mutations and the complete truncation of the carboxy terminus), could be coimmunoprecipitated with WT CFTR. The biochemical studies were complemented by fluorescence lifetime imaging microscopy (FLIM), a distance-dependent approach that monitors protein-protein interactions between fluorescently labeled molecules. Our measurements revealed significantly increased fluorescence resonance energy transfer between CFTR and all tested ENaC combinations as compared with controls (ECFP and EYFP cotransfected cells). Our findings are consistent with the notion that CFTR and ENaC are within reach of each other even in the setting of Liddle's syndrome mutations, suggestive of a direct intermolecular interaction between these two proteins.

Single-cell transcriptomics trajectory and molecular convergence of clinically relevant mutations in Brugada syndrome.

Brugada syndrome (BrS) is a rare, inherited arrhythmia with high risk of sudden cardiac death. To evaluate the molecular convergence of clinically relevant mutations and to identify developmental cardiac cell types that are associated with BrS etiology, we collected 733 mutations represented by 16 sodium, calcium, potassium channels, and regulatory and structural genes related to BrS. Among the clinically relevant mutations, 266 are unique singletons and 88 mutations are recurrent. We observed an over-representation of clinically relevant mutations (∼80%) in SCN5A gene and also identified several candidate genes, including GPD1L, TRPM4, and SCN10A. Furthermore, protein domain enrichment analysis revealed that a large proportion of the mutations impacted ion transport domains in multiple genes, including SCN5A, TRPM4, and SCN10A. A comparative protein domain analysis of SCN5A further established a significant (P = 0.04) enrichment of clinically relevant mutations within ion transport domain, including a significant (P = 0.02) mutation hotspot within 1321-1380 residue. The enrichment of clinically relevant mutations within SCN5A ion transport domain is stronger (P = 0.00003) among early onset of BrS. Our spatiotemporal cellular heart developmental (prenatal to adult) trajectory analysis applying single-cell transcriptome identified the most frequently BrS-mutated genes (SCN5A and GPD1L) are significantly upregulated in the prenatal cardiomyocytes. A more restrictive cellular expression trajectory is prominent in the adult heart ventricular cardiomyocytes compared to prenatal. Our study suggests that genomic and proteomic hotspots in BrS converge into ion transport pathway and cardiomyocyte as a major BrS-associated cell type that provides insight into the complex genetic etiology of BrS.NEW & NOTEWORTHY Brugada syndrome is a rare inherited arrhythmia with high risk of sudden cardiac death. We present the findings for a molecular convergence of clinically relevant mutations and identification of a single-cell transcriptome-derived cardiac cell types that are associated with the etiology of BrS. Our study suggests that genomic and proteomic hotspots in BrS converge into ion transport pathway and cardiomyocyte as a major BrS-associated cell type that provides insight into the complex genetic etiology of BrS.

Long-Read Sequencing Improves the Detection of Structural Variations Impacting Complex Non-Coding Elements of the Genome.

The advent of long-read sequencing offers a new assessment method of detecting genomic structural variation (SV) in numerous rare genetic diseases. For autism spectrum disorders (ASD) cases where pathogenic variants fail to be found in the protein-coding genic regions along chromosomes, we proposed a scalable workflow to characterize the risk factor of SVs impacting non-coding elements of the genome. We applied whole-genome sequencing on an Emirati family having three children with ASD using long and short-read sequencing technology. A series of analytical pipelines were established to identify a set of SVs with high sensitivity and specificity. At 15-fold coverage, we observed that long-read sequencing technology (987 variants) detected a significantly higher number of SVs when compared to variants detected using short-read technology (509 variants) (p-value < 1.1020 × 10-57). Further comparison showed 97.9% of long-read sequencing variants were spanning within the 1-100 kb size range (p-value < 9.080 × 10-67) and impacting over 5000 genes. Moreover, long-read variants detected 604 non-coding RNAs (p-value < 9.02 × 10-9), comprising 58% microRNA, 31.9% lncRNA, and 9.1% snoRNA. Even at low coverage, long-read sequencing has shown to be a reliable technology in detecting SVs impacting complex elements of the genome.

SARS-CoV-2 may hijack GPCR signaling pathways to dysregulate lung ion and fluid transport.

The tropism of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a virus responsible for the ongoing coronavirus disease 2019 (COVID-19) pandemic, toward the host cells is determined, at least in part, by the expression and distribution of its cell surface receptor, angiotensin-converting enzyme 2 (ACE2). The virus further exploits the host cellular machinery to gain access into the cells; its spike protein is cleaved by a host cell surface transmembrane serine protease 2 (TMPRSS2) shortly after binding ACE2, followed by its proteolytic activation at a furin cleavage site. The virus primarily targets the epithelium of the respiratory tract, which is covered by a tightly regulated airway surface liquid (ASL) layer that serves as a primary defense mechanism against respiratory pathogens. The volume and viscosity of this fluid layer is regulated and maintained by a coordinated function of different transport pathways in the respiratory epithelium. We argue that SARS-CoV-2 may potentially alter evolutionary conserved second-messenger signaling cascades via activation of G protein-coupled receptors (GPCRs) or by directly modulating G protein signaling. Such signaling may in turn adversely modulate transepithelial transport processes, especially those involving cystic fibrosis transmembrane conductance regulator (CFTR) and epithelial Na+ channel (ENaC), thereby shifting the delicate balance between anion secretion and sodium absorption, which controls homeostasis of this fluid layer. As a result, activation of the secretory pathways including CFTR-mediated Cl- transport may overwhelm the absorptive pathways, such as ENaC-dependent Na+ uptake, and initiate a pathophysiological cascade leading to lung edema, one of the most serious and potentially deadly clinical manifestations of COVID-19.

Novel mutations in actionable breast cancer genes by targeted sequencing in an ethnically homogenous cohort.

BACKGROUND: Genetic testing is becoming an essential tool for breast cancer (BC) diagnosis and treatment pathway, and particularly important for early detection and cancer prevention. The purpose of this study was to explore the diagnostic yield of targeted sequencing of the high priority BC genes. METHODS: We have utilized a cost-effective targeted sequencing approach of high priority actionable BC genes (BRCA1, BRCA2, ERBB2 and TP53) in a homogeneous patient cohort from Bangladesh (n = 52) by using tumor and blood samples. RESULTS: Blood derived targeted sequencing revealed 25.58% (11/43) clinically relevant mutations (both pathogenic and variants of uncertain significance (VUS)), with 13.95% (6/43) of samples carrying a pathogenic mutations. We have identified and validated five novel pathogenic germline mutations in this cohort, comprising of two frameshift deletions in BRCA2, and missense mutations in BRCA1, BRCA2 and ERBB2 gene respectively. Furthermore, we have identified three pathogenic mutations and a VUS within three tumor samples, including a sample carrying pathogenic mutations impacting both TP53 (c.322dupG; a novel frameshift insertion) and BRCA1 genes (c.116G > A). 22% of tissue samples had a clinically relevant TP53 mutation. Although the cohort is small, we have found pathogenic mutations to be enriched in BRCA2 (9.30%, 4/43) compare to BRCA1 (4.65%, 2/43). The frequency of germline VUS mutations found to be similar in both BRCA1 (4.65%; 2/43) and BRCA2 (4.65%; 2/43) compared to ERBB2 (2.32%; 1/43). CONCLUSIONS: This is the first genetic study of BC predisposition genes in this population, implies that genetic screening through targeted sequencing can detect clinically significant and actionable BC-relevant mutations.

Electrochemical impedance spectroscopy for black lipid membranes fused with channel protein supported on solid-state nanopore.

Black lipid membranes (BLMs) have been used for detecting single-channel activities of pore-forming peptides and ion channels. However, the short lifetimes and poor mechanical stability of suspended bilayers limit their applications in high throughput electrophysiological experiments. In this work, we present a synthetic solid-state nanopore functionalized with BLM fused with channel protein. A nanopore with diameter of ~180 nm was electrochemically fabricated in a thin silicon membrane. Folding and painting techniques were demonstrated for production of stable suspended BLMs followed by incorporation of transmembrane protein, ENaC. Membrane formation was confirmed by employing electrochemical impedance spectroscopy (EIS) in the frequency regime of 10-2-105 Hz. Results show that electrochemically fabricated solid state nanopore support resulted in excellent membrane stability, with >1 GΩ of up to 72 and 41 h for painting and folding techniques, respectively. After fusion of ENaC channel protein, the BLM exhibits the stability of ~5 h. We anticipate that such a solid-state nanopore with diameter in the range of 150-200 nm and thickness <1 µm could be a potential platform to enhance the throughput of ion-channel characterization using BLMs.

Low temperature and chemical rescue affect molecular proximity of DeltaF508-cystic fibrosis transmembrane conductance regulator (CFTR) and epithelial sodium channel (ENaC).

An imbalance of chloride and sodium ion transport in several epithelia is a feature of cystic fibrosis (CF), an inherited disease that is a consequence of mutations in the cftr gene. The cftr gene codes for a Cl(-) channel, the cystic fibrosis transmembrane conductance regulator (CFTR). Some mutations in this gene cause the balance between Cl(-) secretion and Na(+) absorption to be disturbed in the airways; Cl(-) secretion is impaired, whereas Na(+) absorption is elevated. Enhanced Na(+) absorption through the epithelial sodium channel (ENaC) is attributed to the failure of mutated CFTR to restrict ENaC-mediated Na(+) transport. The mechanism of this regulation is controversial. Recently, we have found evidence for a close association of wild type (WT) CFTR and WT ENaC, further underscoring the role of ENaC along with CFTR in the pathophysiology of CF airway disease. In this study, we have examined the association of ENaC subunits with mutated ΔF508-CFTR, the most common mutation in CF. Deletion of phenylalanine at position 508 (ΔF508) prevents proper processing and targeting of CFTR to the plasma membrane. When ΔF508-CFTR and ENaC subunits were co-expressed in HEK293T cells, we found that individual ENaC subunits could be co-immunoprecipitated with ΔF508-CFTR, much like WT CFTR. However, when we evaluated the ΔF508-CFTR and ENaC association using fluorescence resonance energy transfer (FRET), FRET efficiencies were not significantly different from negative controls, suggesting that ΔF508-CFTR and ENaC are not in close proximity to each other under basal conditions. However, with partial correction of ΔF508-CFTR misprocessing by low temperature and chemical rescue, leading to surface expression as assessed by total internal reflection fluorescence (TIRF) microscopy, we observed a positive FRET signal. Our findings suggest that the ΔF508 mutation alters the close association of CFTR and ENaC.

CFTR regulation of epithelial sodium channel.

Cystic fibrosis (CF) is a lethal genetic disorder, characterized by both clinical and genetic complexities, and arises as a result of mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene. The gene encodes a Cl(-) channel belonging to the ABC (ATP Binding Cassette) superfamily of transporters. The members of this superfamily use ATP hydrolysis to fulfill their function as active transporters. So far, CFTR is the only member of this family to function as a cAMP-activated Cl(-) channel. Intense research following the cloning of the CFTR gene has extended the role of the CFTR beyond that of a Cl(-) channel. One of the best recognized, yet still controversial, functions of the CFTR is its ability to modulate the functioning of other transporters. The modulation of epithelial Na(+) channel (ENaC) function serves as a prime example of regulatory function of the CFTR. In this chapter, we will briefly describe an integrated protocol consisting of biochemical and electrophysiological approaches to study the regulation of ENaC by CFTR.

Psalmotoxin-1 docking to human acid-sensing ion channel-1.

Acid-sensing ion channel-1 (ASIC-1) is a proton-gated ion channel implicated in nociception and neuronal death during ischemia. Recently the first crystal structure of a chicken ASIC was obtained. Expanding upon this work, homology models of the human ASICs were constructed and evaluated. Energy-minimized structures were tested for validity by in silico docking of the models to psalmotoxin-1, which potently inhibits ASIC-1 and not other members of the family. The data are consistent with prior radioligand binding and functional assays while also explaining the selectivity of PcTX-1 for homomeric hASIC-1a. Binding energy calculations suggest that the toxin and channel create a complex that is more stable than the channel alone. The binding is dominated by the coulombic contributions, which account for why the toxin-channel interaction is not observed at low pH. The computational data were experimentally verified with single channel and whole-cell electrophysiological studies. These validated models should allow for the rational design of specific and potent peptidomimetic compounds that may be useful for the treatment of pain or ischemic stroke.

Assessment of the CFTR and ENaC association.

Cystic fibrosis (CF) is one of the most common lethal genetic disorders. It results primarily from mutations in the cystic fibrosis transmembrane conductance regulator (cftr) gene. These mutations cause inadequate functioning of CFTR, which in turn leads to the severe disruption of transport function in several epithelia across various organs. Affected organs include the sweat glands, the intestine, and the reproductive system, with the most devastating consequences due to the effects of the disease on airways. Despite aggressive treatment, gradual lung failure is the major life limiting factor in patients with CF. Understanding of the exact manner by which defects in the CFTR lead to lung failure is thus critical. In the CF airway, decreased chloride secretion and increased salt absorption is observed. The decreased chloride secretion appears to be a direct consequence of defective CFTR; however, the increased salt absorption is believed to result from the failure of CFTR to restrict salt absorption through a sodium channel named the epithelial Na(+) channel, ENaC. The mechanism by which CFTR modulates the function of ENaC proteins is still obscure and somewhat controversial. In this short review we will focus on recent findings of a possible direct CFTR and ENaC association.

Molecular proximity of cystic fibrosis transmembrane conductance regulator and epithelial sodium channel assessed by fluorescence resonance energy transfer.

We present the evidence for a direct physical association of cystic fibrosis transmembrane conductance regulator (CFTR) and epithelial sodium channel (ENaC), two major ion channels implicated in the pathophysiology of cystic fibrosis, a devastating inherited disease. We employed fluorescence resonance energy transfer, a distance-dependent imaging technique with capability to detect molecular complexes with near angstrom resolution, to estimate the proximity of CFTR and ENaC, an essential variable for possible physical interaction to occur. Fluorescence resonance energy transfer studies were complemented with a classic biochemical approach: coimmunoprecipitation. Our results place CFTR and ENaC within reach of each other, suggestive of a direct interaction between these two proteins.

Epithelial sodium channel in planar lipid bilayers.

Amiloride-sensitive Na+ channels belong to the epithelial Na+ channel (ENaC)-degenerin superfamily of ion channels. In addition to their key role in sodium handling, they serve diverse functions in many tissues. Improper functioning of ENaC has been implicated in several diseases, including salt-sensitive hypertension (Liddle's syndrome), salt-wasting syndrome (pseudohypoaldosteronism type I), pulmonary edema, and cystic fibrosis. We have utilized planar lipid bilayers, a well-defined system that allows simultaneous control of "internal" and "external" solutions, to study ENaCs.

Methods to study CFTR protein in vitro.

CFTR is a cyclic AMP and nucleotide-related chloride-selective channel with a low unitary conductance. Many of the physiological roles of CFTR are effectively studied in intact cells and tissues. However, there are also several clear advantages to the application of cell-free technologies to the study of the biochemical and biophysical properties of CFTR. When expressed in heterologous cells, CFTR is processed relatively poorly, depending, however, on the cell-type analysed. In some cells, only 20-25% of the protein which is initially synthesized exits the endoplasmic reticulum to insert into the cell membrane [Cell 83 (1995) 121; EMBO J. 13 (1994) 6076]. Further, many of the disease-causing mutants of CFTR result in even lower processing efficiencies. Therefore, several procedures have been developed to study regulated CFTR channel function expressed in microsomal membranes and following its purification and reconstitution. These experimental approaches and their application are discussed here.